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ABSTRACT 



k> ' Aims. We present early optical spectroscopy of the afterglow of the gamma-ray burst GRB 060206 with the aim of determining the metallicity 
of the GRB absorber and the physical conditions in the circumburst medium. We also discuss how GRBs may be important complementary 
probes of cosmic chemical evolution. 

Methods. Absorption line study of the GRB afterglow spectrum. 

Results. We determine the redshift of the GRB to be z = 4.04795 ± 0.00020. Based on the measurement of the neutral hydrogen column 
density from the damped Lyman-o- line and the metal content from weak, unsaturated S ii lines we derive a metallicity of [S/H]= -0.84 ± 0. 10. 
This is one of the highest metallicities measured from absorption lines at z ~ 4. From the very high column densities for the forbidden Si ii*, 
Oi*, and Oi** lines we infer very high densities and low temperatures in the system. There is evidence for the presence of H2 molecules 
with log yV(H2)~17.0, translating into a molecular fraction of log/ a -3.5 with / = 2A'(H2)/(2A'(H2)+ A'(Hl)). Even if GRBs are only 
formed by single massive stars with metallicities below ~ 0.3 Zq, they could still be fairly unbiased tracers of the bulk of the star formation at 
z > 2. Hence, metallicities as derived for GRB 060206 here for a complete sample of GRB afterglows will directly show the distribution of 
metallicities for representative star-forming galaxies at these redshifts. 



Key words. Gamma rays: bursts - galaxies: high redshift, abundances - cosmology: observations 



1. Introduction 

Long gamma-ray bursts (GRBs) are now established to be 
caused by the deaths of massive stars (e.g., ,Hjorth et aL 2003. 
IStanek et aL 20031 and due to their brightness they can be 



* Based on observations made at the Spanish Observatorio del 
Roque de los Muchachos on La Palma with the Nordic Optical 
Telescope, operated jointly by Denmark, Finland, Iceland, Norway, 
and Sweden, and with the William Herschel Telescope. 



observed throughout most of the observable Universe (e.g., 
'Kawai et a L 2006t . Given these facts it has long been real- 
ized that GRBs could be ideal probes of star-formation activity 
throughout the history of the Universe (e.g., |Wijers et al. 1998^ . 
However, previous GRB missions detected too few rapidly 
well-localised GRBs to build statistically interesting samples 
and hence really capitalize on this potential. The Swift satellite 
(Gehrels et al. 2004 1 has increased the detection rate of rapidly 
well-localised GRBs by roughly an order of magnitude com- 
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pared to previous missions. Moreover, its significantly deeper 
detection limit (e.g., Band i2006t means that Swift detects more 
distant bursts than previous missions f Jakobsso n et al. 2006> . 

In this Letter we present optical spectroscopy of 
GRB 060206, focussing on the measurement of the metallicity 
of the GRB absorption system, the density in the circumburst 
environment and the presence of molecules. We then discuss 
how GRBs are important for understanding the build-up and 
distribution of metals in galaxies (e.g., Pei & Fallll995^, which 
is currently not fully understood (e.g., Ferrara et al. .2005i . 

2. Observations 

GRB 060206 was discovered by the Burst Alert Telescope 
(BAT) aboard the Swift satelUte on February 6 04:46:53 UT. 
The burst exhibited a slow rise and a faster decline, with a 
Tgo of 7 + 2 s (Palmer et al. I2006> . The X-ray Telescope 
(XRT) slewed promptly to the location and began taking data 
at At = 58 s, where Af is the time from the onset of the burst. 
Due to entry into the South Atlantic Anomaly, XRT could only 
observe the BAT error circle briefly, and therefore no fading 
X-ray source was immediately localised (Morris et al. 2006a i. 
Observations with the UV/Optical Telescope (UVOT) began at 
Af = 57 s but failed to reveal an optical afterglow (OA) candi- 
date in the initial data products (Morris et al. 2006a i. 

We observed the GRB 060206 BAT error circle in the 
R-bwd with the Andalucia Faint Object Spectrograph and 
Camera (ALFOSC) on the Nordic Optical Telescope (NOT) 
starting at Af » 15min. A point-like object (R ~ 17.3) not 
present in the Digitized Sky Survey was detected. The detec- 
tion was confirmed by re-analysis of the XRT and UVOT data 
(Boyd et al. 2006, Morris et al. 2006b i. 

Starting at Af x 48 min we obtained a 1800 s spectrum with 
a low resolution (LR) grism and a 1.3 arcsec wide slit covering 
the spectral range from about 3500 A to 9000 A at a resolution 
of 14 A. The airmass during the spectroscopic observation was 
very low, starting at 1.01. We also obtained medium resolu- 
tion (MR) spectroscopic data with the Intermediate-dispersion 
Spectroscopic and Imaging System (ISIS) on the William 
Herschel Telescope (WHT) at Af = 1.61 hours. We took spec- 
tra of exposure times 2x900 s with both the the 600B grating 
(covering 3800-5300 A) and and the 1200R grating (covering 
6100-7200 A) both with a 1.0 arcsec wide slit. The mean air- 
mass was 1.036 and the observations were carried out at the 
parallactic angle. The resolutions of blue- and red-arm spectra 
are 1.68 and 0.82 A, respectively. 

3. Results 

The spectra show absorption both at the GRB redshift of z = 
4.04795 + 0.00020 and from two intervening systems at red- 
shifts z = 1.48 and z = 2.26. In the analysis presented here 
we focus on the metallicity of the GRB absorption system. 
From a fit of the damped Lya line at z = 4.048, we measure a 
neutral hydrogen column density of log A^(Hi) = 20.85 ± 0.10 
(Fig. [3, well above the definition for Damped Lya Absorbers 
(DLAs, logA^(Hi) > 20.3). This value is consistent with fits 
to the higher order Lyman lines in the blue MR spectrum. To 
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Fig. 1. Upper panel: Section of the LR afterglow spectrum 
showing the DLA line at the GRB redshift, Zabs - 4.048. 
Overlaid is the best fitting DLA profile, corresponding to 
logA^(Hi) = 20.85 + 0.10. Middle panel: Fits to the Oi, Oi*, 
Oi**, Sin, Sin*, and S ii lines in the MR spectrum. The zero- 
point for the velocity scale is Zabs = 4.048. Lower panel: The 
likely H2 lines at the GRB redshift in the MR spectrum. 



derive the metal column density we use the unsaturated S n 
/11250, 1253, 1259 lines shown in Fig. □ The profile of the 
GRB absorption systems consists of at least four components 
spread over ~500 km s ' in velocity space. On the right hand 
side of the middle panel of Fig.^ a four-component fit of S II 
and Si II* is shown. On the left hand side, we have only fitted 
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Table 1. Ionic column densities in individual components of 
the GRB system at z^bs - 4.048. 



Ion 


Transition 


log A' ± CTiogA, 


& ± cTfc 




lines used 




(km s ) 


Zabs = 4.0441 










Sii 


1250, 


1253, 1259 


14.36±0.06 


14±7 


Si 11*^ 


1309 




< ij.lU 


t4± / 


Zab.s = 4.0480 










Sii 


1250, 


1253, 1259 


14.79±0.05 


~ 15 


Sin* 


1309 




13.54±0.08 


~ 15 


Zabs = 4.0490 










Sii 


1250, 


1253, 1259 


14.60±0.06 


~ 15 


Sin* 


1309 




13.77±0.06 


~ 15 


Zabs = 4.0509 










Sn 


1250, 


1253, 1259 


14.59±0.06 


14±2 


Sin* 


1309 




14.22±0.05 


14±2 


Oi* 


1304 




15.03±0.16 


14±2 


1** 


1306 




14.25±0.05 


14±2 


Sin 


1304 




14.56±0.10 


14±2 


" 3cr upper limit. 



the redmost component at 175 km & where we can derive a 
column density for Si II, O I* and O I**. We fit the lines using 
the FitLyman package in MIDAS; the measured column den- 
sities for all components are given in Tabled Note that for 
components 2 and 3 we fixed the turbulent broadening param- 
eter value to 15 km s"' . Summing over the four components of 
the system, we get logA^(SII) = 15.21 + 0.03 which leads to 
[S/H] = -0.84 + 0.10. 

We also find logA^(SiII*) = 14.42 ± 0.02 and therefore 
[SiII*/SiII]= -1.15, assuming [S/Si] = 0. In the fourth com- 
ponent, we have [SiII*/SiII] - -0.34, which is high com- 
pared to the ratios found in other GRB host galaxies (twice as 
high as the value found for GRB 050505 by Berger et al. l2005l 
and 20 times higher than the value found for GRB 030323 by 
Vreeswijk et al. 2004 1. Fine-structure levels can be populated 
through collisions, photo-excitation by IR photons, and/or flu- 
orescence (Bahcall & Wolf 1968). Assuming that the former 
mechanism is dominant (but see Berger et al. 2006; Prochaska 
et al. 2006), we can estimate the H I volume density, using the 
calculations by Silva & Viegas (2002). Assuming an electron 
fraction He ~ 10"'*«hi (see Vreeswijk et al. 2004; Berger et 
al. 2005) and a temperature of 1000 K, we find «hi ~ 10^ cm"-' 
for the fourth component, and riui ~ 10* cm"^ for the mean 
ratio of all components (see Fig. 8 of Silva & Viegas 2002). 
For the fourth component, we can actually constrain the ki- 
netic temperature and volume density, as the observed ratio 
[OP/OI**] = 0.78 + 0.17 can only be reached below 320 K 
(wifliin the Icr errors), while [SiII*/SiII] = -0.34 ± 0.11 re- 
quires a temperature above 240 K. This temperature range cor- 
responds to a volume density of net = 1-3x10^ cm"-'. 

Finally, our data show the first evidence for H2 molecules 
in a GRB absorber. Two consistent features are seen at the lo- 
cation of the Wl-0 R(0), Wl-0 R(l) and Wl-0 Q(l) lines at 
Zabs - 4.04793 (bottom panel in Fig. with column densi- 
ties of log A^(H2) ~ 17.0 for the J = 1 rotational level, and 



Table 2. Table of pubUshed absorption metallicities for GRBs. 
References: [1] Savaglio et al. (2003); [2] Vreeswijk et 
al. (2006); [3] Vreeswijk et al. (2004); [4] Watson et al. (2006); 
[5] Berger et al. (2005); [6] Starling et al. (2005); [7] Chen et 
al. (2005); [8] Ledoux et al., 2005; [9] Kawai et al. (2005); [10] 
this work. 



GRB 


Metallicity 


Redshift 


Ref 


000926 


[Zn/H]= -0.13 ±0.25 


2.038 


1 


011211 


[Si/H]= -0.90 + 0.5 


2.142 


2 


030323 


[S/H]= -1.26 ±0.2 


3.372 


3 


050401 


[Zn/H]= -1.0 ±0.4 


2.899 


4 


050505 


[S/H]> -1.2 


4.275 


5 


050730 


[S/H]= -2.0 ±0.2 


3.968 


6,7 


050820 


[Si/H]= -0.6 ±0.1 


2.615 


8 


050904 


[S/H]= -1.3 ±0.3 


6.295 


9 


060206 


[S/H]= -0.84 ±0.10 


4.048 


10 



log A^(H2) < 16.7 for J = 0. The corresponding H2 molecular 
fraction is log/ ~ -3.5 with / = 2N(H2)/(2N(¥l2)+ N(¥li). 
This is the second highest redshift at which these H2 lines have 
been detected (Ledoux et aI. I2006> . 

4. Discussion and conclusions 

In Table |5] we have compiled metallicity measurements for 
z > 2 GRB absorption systems from the literature. As seen, 
despite its high redshift, the GRB 060206 system has one of 
the highest metallicities measured for a GRB absorption sys- 
tem and one of the highest metallicities measured from QSO 
absorption lines at z > 4. We also note that the metallicity of 
GRB 060206 is about 15 times higher than that of GRB 050730 
which has a similar redshift (z = 3.968). This shows that GRBs 
do not only occur in very low metallicity environments, but 
also in environments covering a broad range of metallicities at 
a given redshift. 

Fine-structure fines are ubiquitous in GRB absorbers 
(Vreeswijk et aL EUO?! Berger et al. 2005 Chen et al. ,2005t . 
As discussed above, for GRB 060206 we derive high densities 
and low temperatures, consistent with the possible detection of 

H2 in the spectrum, with a molecular fraction of log/ 3.5. 

Molecules are also detected in QSO-DLAs (Ledoux et aI. I2003> 
in cold gas with T< 300 K, but with lower densities (< 400 
cm"^, Srianand et al. .2005.) . The difference in densities along 
random (cross-section selected) and GRB sightlines like the 
one towards GRB 060206 must reflect the higher than aver- 
age density of the star forming region in which the progeni- 
tor star was located. Such regions must have a cross-section 
which is less than a few percent of the total cross-section for 
QSO-DLAs at redshifts z - 2-4 in order to explain that similar 
fine-structure lines have not yet been seen in QSO-DLAs. 

In the collapsar models GRBs can only be formed by mas- 
sive single stars with a metallicity below ~0.3 Zq (Hirschi et 
al. 120051 Woosley & Heger 2006). Such a bias is not present 
in all models, most noteably not in models involving a binary 
progenitor (Fryer & Heger 2005 1. To gauge what a metallicity 
bias could mean for the completeness of GRBs as cosmolog- 
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ical probes it is important to know the present mass-weighted 
mean metallicity and how it declines with redshift. Zwaan et 
al. J2005 I) find a present-day mean metallicity in the gas phase 
of Z a; 0.44 Z0 and a slope between -0.25 and -0.30 dex per 
unit redshift. This means that at z ~ I and earlier the mean 
metallicity of the gas is below the cut-off value above which 
single massive stars in the coUapsar models do not make GRBs. 
It is therefore likely that GRBs at z > 2 will be fairly unbiased 
tracers of star-formation, while they become increasingly bi- 
ased at z < 1 if there is a low-metallicity bias. 

Metallicities can be measured from optical spectroscopy of 
the afterglow or the host galaxy. Fig.|2shows that most metal- 
licities for ORB absorption systems fall below 0.3 Z©, but one 
(GRB 000926, Z = 0.7 Zq, Savagho et al. l2003t has a metal- 
licity well above the threshold for a single massive star to pro- 
duce GRBs in the collapsar models. The metallicity of the host 
galaxy of GRB 980425 by Sollerman et al. ( 2005} is also higher 
than this limit (0.8 Zq). This indicates that collapsars result- 
ing from single massive stars are not the only progenitors to 
long GRBs or that massive stars with Z > 0.3 Zq can also pro- 
duce long GRBs. The GRB metallicities fall within the range 
spanned by QSO-DLAs, but most are above the mean curve 
derived by Zwaan et al. (2005 1. The mean offset is 0.49 and 
the scatter 0.38 dex (excluding GRB 050904 which is at very 
high redshift where there are no QSO-DLA data). This off- 
set may reflect that DLAs are cross-section selected, whereas 
GRBs represent sightlines to more central regions in their hosts 
(Bloom et al. 2002 1. A metallicity gradient of -0.09 dex kpc"' 
as assumed in the study by Zwaan et al. (I2005t (and observed in 
the Galaxy) seems sufficient to explain the offset. Note that the 
current sample of QSO-DLAs is not likely to be strongly biased 
against high metallicity sightlines (Akerman et al. l2005> . 

Furthermore, in the derivation of the cosmic star-formation 
density, and hence the production of metals, the shape of the 
luminosity function is important. Roughly, two thirds of the 
UV light from LBGs is emitted by galaxies that are fainter 
than the actual flux limit of the ground based LBG survey of 
Steidel et al. (2003). Using only LBGs to derive the total star- 
formation density means that an extrapolation to the poorly de- 
termined faint end of the luminosity function is unavoidable. 
GRBs allow us to probe this faint end since the selection is not 
limited by the brightness of the host. We refer to Jakobsson 
et al. (1200 5 1 for a quantitative analysis. Finally, GRBs allow 
the measurement of metallicities at very high redshifts (z > 6, 
IKawai et al. 2006( l. which is inaccessible to QSO-DLAs. 

In conclusion, we measure a metallicity of [S/H]= -0.84 + 
0.10 for GRB 060206 at z = 4.048. We also find the first ev- 
idence for molecular lines from a GRB absorber We have ar- 
gued that GRBs can be used to measure the metallicities and 
luminosities of typical star forming galaxies at z > 2, making 
GRBs promising complementary probes of chemical evolution 
at high redshift. 
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